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Abstract 

Optical and electrical characteristics of AlGaAs lasers with separate confinement 
heterostructures are modeled by using Synopsys's Sentaurus TCAD, and open source 
software. The results for cases of 2-QW (2 Quantum Wells) and 3-QW structures are 
compared with these for 1-QW. A significant improvement of useful laser parameters 
is obtained by increasing the number of Quantum Wells and optimizing the width 
of waveguides. In particular, the maximum optical efficiency is shown to reach 88 
% for a 3-QW structure with optimal width of waveguides. The width of optical 
intensity profile of MQW lasers increases, leading to lowering maximal light power 
density passing through laser facets, decreasing the risk of catastrophic damage of 
mirrors. 



1 Introduction 

Alferov [JJ, et al., proposed creating semiconductor-based lasers comprising the use 
of a geometrically-narrow active recombination region where photon generation oc- 
curs, with waveguides around improving the gain to loss ratio (separate confinement 
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heterostructures; SCH). That idea dominated largely the field of optoelectronics de- 
velopment in the past years. Due to the relative simplicity and perfection of technol- 
ogy, solid solutions of Al x Ga\- x As are commonly used as wide-gap semiconductors 
in SCH lasers. 

Finding however the best, most efficient structures of SCH lasers is not a task 
that could be solved easy by experimenting only. Computational modeling is being 
used recently, as the cost effective approach towards designing new devices as well 
predicting their characteristics. Moreover, computational modeling offers also new 
approaches to study physical properties of a device: often, limits of technology do 
not allow to carry out experiments with accuracy high enough that would allow to 
verify certain hypothesis. 

In our earlier works we first were able to find agreement between our calculations 
of quantum well energy states and the lasing wavelength observed experimentally [2]. 
Next p], several changes in structure and doping of SCH AlGaAs lasers have been 
shown to considerably improve their electrical and optical parameters. We compared 
computed properties with these of lasers produced by Polyus research institute in 
Moscow [1], [3]. In particular, by changing the width of active region (Quantum 
Well), waveguide width, doping concentration in all laser layers, as well by changing 
the waveguide profile by introducing a gradual change of Al concentration, and also 
by introducing variable doping profile of carriers across waveguide, we were able to 
decrease significantly the lasing threshold current, increase the slope of optical power 
versus current, and increase optical efficiency up to about 74 % [6]. We have shown 
[7] that the lasing action may not occur at certain widths or depths of Quantum 
Well (QW), and the threshold current as a function of these parameters may have 
discontinuities that occur when the most upper quantum well energy values are very 
close to either conduction band or valence band energy offsets. These effects are more 
pronounced at low temperatures, and may be observed also, at certain conditions, 
in temperature dependence of lasing threshold current as well. 

A simple analytical, phenomenological model was shown to describe optical effi- 
ciency, rj = L/P (L is optical power intensity and P is the supplied electrical power), 
with a high accuracy, by using two parameters only [8j. 

In Multi-Quantum Well (MQW) structures better characteristics are achievable 
than in case when one QW is used. Hence, a question is to what an extend yet could 
we improve properties of laser studied, by using optimized results obtained already 
and by introducing MQW. The purpose of this work is to compare results for 1, 2, 
and 3 QW structures and find out optimal widths of waveguides in each case. 

One of the fundamental laser characteristics are their threshold current, Ifh, 
above which lasing action starts, the slope of optical power versus current, S = 
dL/dl, which is approximately constant for currents just above J^, and lasing off- 
set voltage Uq. We choose these parameters as characteristics to be compared in 
modeling. 
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Table 1: Reference structure of AlGaAs 1-QW SCH laser layers used in computer model- 
ing, d is the width of layers. 



No Layer Composition Doping [cm 3 ] d [fim] 



1 


n-substrate 


n-GaAs (100) 


2 


•10 18 


350 


2 


n-buffer 


n-GaAs 


1 


•10 18 


0.4 


3 


n-emitter 


Al _ 5 Ga _ 5 As 


1 


•10 18 


1.6 


4 


waveguide 


Al . 33 Ga . 67 As 


n 


« 10 15 


0.2 


5 


QW 


AlomGaomAs 


n 


« 10 15 


0.012 


6 


waveguide 


Alo.33Gao.e7 As 


n 


« 10 15 


0.2 


7 


p-emitter 


Alo. 5 Ga . 5 As 


1 


•10 18 


1.6 


8 


contact 


p-GaAs 


4 


•10 19 


0.5 



2 Lasers structure and calibration of modeling. 

The reference laser we model has a 1000/xm cavity length and 100/xm width, with 
doping and Al-content as described in Table [TJ The Table [2] describes it's experimen- 
tal parameters. The structure of 2-QW and 3-QW lasers, as shown schematically in 
Figure [H is described in Tables [3] and [4} respectively. 

Synopsys's Sentaurus TCAD is used for modeling [S]. This is an advanced, 
flexible commercial computational environment used for modeling a broad range of 
technological and physical processes in microelectronics world. In case of lasers, some 
calculations in Sentaurus have purely phenomenological nature. The electrical and 
optical characteristics depend, primarily, on the following computational parameters 
that are available for adjusting: 

AreaFactor of electrodes, A e , AreaFactor in Physics section, A p h, electrical con- 
tact resistance R x , and reflection coefficient of laser mirrors, Ri and R r . There are 
several parameters for adjustment that are related to microscopic physical properties 
of materials or structures studied. However, often their values are either unknown 
exactly or finding them would require quantum-mechanical modeling based on first- 
principles. This is however not the aim of our work. 

In order to find agreement between the calculated results and these observed 
experimentally we adjust accordingly values of A e and A p h- 

The results for Ith and S are all shown normalized by the values for the reference 
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Table 2: Summary of experimental conditions and laser parameters, for the reference 



Temperature [K] 


300 


Lasing wavelength [nm] 


808 


Offset voltage U [V] 


1.56-1.60 


Differential resistance, r = dll/dl [mO] 


50-80 


Threshold current I® h [mA] 


200-300 


Slope of optical power, So = dL/dl [W/A] 


1.15-1.25 


Left mirror reflection coefficient Ri 


0.05 


Right mirror reflection coefficient R r 


0.95 



laser described in Table [H 1?. and So, respectively. 

We neglect here the effect of contact resistance, by not including buffer and 
substrate layers and contacts into calculations (compare with structure described in 
Table [TJ. We use instead Inner Voltage parameter available in Sentaurus and treat it 
as a physical quantity that is related to voltage applied. However, we examined in 
details results of calculations of optical and electrical characteristics and compared 
them with these from measurements, finding a reasonably good agreement between 
them, for lasing offset voltage, threshold current, optical intensity, and differential 
resistivity [8]. 

3 Threshold current and dL/dl for 2-QW and 
3-QW structures. 

We choose to use threshold current Ith an d the slope of optical intensity versus 
current, S = dL/dl as the most important technologically laser parameters, when 
comparing lasers with various structures. The advantage of using them is also in 
that that these could be easy extracted semi-automatically from a large collection 
of data sets. Details of data analysis are described better in [8J and also on our 
laboratory web sit^]. 

Computation was performed as a function of the width d w i, for a certain set of 



1 Web: www.ostu.ru/units/ltd/zbigniew/synopsys.php 



4 



Table 3: Two Quantum Wells: Structure of AlGaAs SCH laser layers, d is the width of 
layers. 



AT 

No 


Layer 


Composition 


Dopin 


g [cm \ 


j r 1 
d [/imj 


1 


n-emitter 


Al . 5 Gao.5As 


N 


1 ■ 


10 18 


1.6 


2 


waveguide 


Alo.33Gao.67As 


N 


1 • 


10 15 


d w o 


3 


QW 


Al0.0sGa0.g2 As 


N 


1 ■ 


10 15 


0.012 


4 


waveguide 


Alo.33Gao.e7 As 


N 


5- 


10 14 


d w i 


5 


QW 


AlomGao. 92 As 


P 


1 • 


10 15 


0.012 


6 


waveguide 


Alo.33Gao.67 As 


P 


1 • 


10 15 


d w o 


7 


p-emitter 


Alo. 5 Ga . 5 As 


P 


1 • 


10 18 


1.6 



constant values of d w o (see Fig. [I]). 

Figures [2] and [H] show results for J^/J^ and S/So, respectively, for the case of 
2-QW structures, while Figures U] and [5] show similar results for 3-QW structures. 
In these Figures, where no data points are present for certain values of d w \ - there is 
no lasing action observed there. The meaning of lines in Figures [2][1] is to guide the 
eyes only. They are drawn however by using least-squares fitting of the data points 
to simple power-law functions. 

As we see, a significant reduction of threshold current is observed for both, 2-QW 
and 3-QW structures, if compared with the value for the reference laser with 1-QW. 
Also, a large increase of dL/dl is found, more significant in case of 3-QW structure. 

4 Optical efficiency. 

Figures [6] and [7J illustrate how optical efficiency as a function of voltage applied 
depends on geometry of lasers, for several cases of waveguide parameters d w o and 
d w \ that correspond, approximately, to the maximum values of optical efficiency. 
The largest optical efficiency achieved, of about 88%, is for 3-QW structure when 
d w o is about 100 nm, and d w \ of about 29 nm. 

We proposed [8J that a modified exponential dependence describes very well 
current-voltage characteristics just above the lasing offset U > Uo'- 

I(U) = I th ■ exp(C ■ (U - Uo) + D ■ (U - U ) 2 ), (1) 
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Table 4: Three Quantum Wells: Structure of AlGaAs SCH laser layers, d is the width 
layers. 



No Layer Composition Doping [cm 3 ] d [fim] 



1 


n-emitter 


Al . 5 Gao.5As 


N 


1 


■10 18 


1.6 


2 


waveguide 


Air, wGdn 67 As 


N 


1 


•10 15 


d<wQ 


3 


QW 


Al0.0sGa0.g2 As 


N 


1 


■10 15 


0.012 


4 


waveguide 


Alo.33Gao.e7 As 


N 


5 


■10 14 


d w i 


5 


QW 


AlomGao. 92 As 


N 


5 


•10 14 


0.012 


6 


waveguide 


Alo.33Gao.67 As 


P 


5 ■ 


10 14 


d w i 


7 


QW 


Alo.08Gao.g2As 


P 


5- 


10 14 


0.012 


8 


waveguide 


Alo.33Gao.67 As 


P 


1 • 


10 15 


d w o 


9 


p-emitter 


Alo.5Gao.5As 


P 


1 • 


10 18 


1.6 



where C and D are certain fitting parameters. It is convenient to rewrite [T] in 
dimensionless variables: 



i{u) = exp (4 • (u - 1) • [1 + P ■ f/ 2 ■(«-!)]), 



(2) 



where we defined: i(u) = I{U)/Ith and u = U/Uq, a = r ■ Ith/Uo, /3 = Uo • D/C, 
and we defined also r = dU/dl = I 1 g , which is differential resistivity just above 

Hence, the optical efficiency, rj = L/(U ■ I), is given by: 



n ( u \ - JL . M=l (3) 

''W U u-i(u) ' 

where i(u) is given by Eq. [2J 

Like in case of 1-QW structure ([8j), the data shown in Figures [6] and [7] are 
described by [3] with a high accuracy, by using two fitting parameters only, a and f3 
(except of Ith and Uq that may be found in a straightforward way). 
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n-emitter 



d w o d wl d w0 



p-emitter 



d w o d wl d wl d w0 



Figure 1: Schematic structure of energy gap for MQW lasers modeled, with 2- and 3- 
quantum wells (the upper and lower diagrams, respectively). d W Q and d w \ are widths of 
waveguide regions. Active regions remain the same and their width is Ylnm. Doping 
types and concentrations are shown in Tables |3] and HI 



5 Lasing offset voltage U { 



CI- 
As shown by Eq. [3j optical efficiency depends on lasing offset voltage and it is 
desirable to have it's value as low as possible. In literature, there is no a simple 
formula that would describe variations of Uq with physical material parameters and 
geometrical structure of a laser. One of reasons for that is that the I — V charac- 
teristics depend strongly in a too complex way on density of states and number of 
captured carriers in QW. The number of QW levels and their separation, as well 
their energy distance from conduction band or valance band offsets are a function of 
QW geometry as well [7J. Therefore, it is useful to have an insight into what kind of 
dependencies might be obtained as a function of waveguide width. Performing real 
experiments of that kind is unrealistic: there is no way to control laser geometrical 
dimensions during technological process with sufficient, required accuracy. Instead, 
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0.2 - 



1 1 1 1 1 

0.05 0.1 0.15 0.2 

dwj [Mm] 

Figure 2: Two Quantum Wells: threshold current, I t h, normalized by I® h , as a function of 
the width of waveguide in the center, for a few width values of other two waveguides (in 
nm), as indicated in the Figure. 



a large spread of the data points would be obtained and no conclusions could be 
drawrH. 

We find in our modeling that a good accuracy of determining Uq is gain versus 
voltage curve in a near range of voltage values below Uq. We use a linear extrap- 
olation of the data towards the maximum value, which is constant above Uq in 
Sentaurus. The results presented in Fig. [8] were obtained that way. This method is 
the most accurate in our case and allows for easy, semi-automatic analysis of large 
collections of data sets. 

As a result of a complex interplay between density of states in quantum wells 
when waveguide width d w \ changes, clear discontinuities are observed in Uo(d w i), 
as shown on Fig. [SI 

2 In [7] , we point out however that there are hypothetically ways that could bypass these experimental 
problems. 



8 




1.05 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 

d w i [Mm] 

Figure 3: Two Quantum Wells: Slope of optical power, S = dL/dl, normalized by S , 
corresponding to the data in Figure [2J 



One should however be cautioned here. Somewhat similar effects might be an 
undesired, artificial artifact of improperly conducting the modeling calculations: 
results depend for instance on how the mesh changes when geometrical parameters of 
modeled structure change. In our case, we kept the number of mesh points constant, 
independent of geometrical sizes of laser layers. 



6 Optical intensity profiles for 3-QW structure. 

When changing the structure and width of waveguides we gain some control over 
the power distribution over the laser facets. That way we may possibly decrease the 
effects that limit the maximal output power, that are caused by thermal damage to 
mirrors. 

One of the simplest ways to characterize optical intensity profiles is to use their 
half-width, i.e. the distance from the laser center where optical intensity decreases 
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0.6 1 1 1 1 1 

0.05 0.1 0.15 0.2 

d w i [nm] 



Figure 4: Three Quantum Wells: threshold current, I t h, normalized by I® h , as a function 
of the width of two waveguides in the center, for a few values of other two waveguides 
width (in nm), as indicated in the Figure. 



to half of its maximal value, Aw. 

The half-width in Figure [9J is given by linear dependence on Wq: Aw = 0.409 • 
u>o + 0.161 (when w\ = 0.130//m). 

When similar analysis is done for the case of w\ width of 0.029fim, we obtain 
the following expression for the half-width of optical profile intensity: Aw = 0.394 • 
w + 0.111. 

For the case of w\ width of 0.219^m, we obtain the following expression: Aw = 
0.444 • w + 0.2034. 

By additional analysis, we find that the following general relation describes well 
the dependence of the half- width of optical intensity profile on wo and w\ : 

Aw = (0.392 + 1.069 • w\) ■ w + 0.487 • w x + 0.097 (4) 
where all quantities are in \xm. 
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c 
55 




0.05 0.1 0.15 0.2 

d w i [^ m ] 

Figure 5: Three Quantum Wells: Slope of optical power, S = dL/dl, normalized by So, 
corresponding to the data in Figure HI 

7 Summary and Conclusions. 

By using Synopsys's Sentaurus TCAD, and open source software we performed com- 
puter modeling of optical and electrical characteristics of AlGaAs lasers with sep- 
arate confinement heterostructures, when 2 and 3 quantum wells are present. We 
compared results with these for 1-QW laser calibrated to reproduce characteristics 
of lasers produced at Polyus research institute in Moscow ([1] and [5j). 

It was shown that a significant improvement of laser parameters is obtained in 
case of MQW lasers, when the width of waveguides is optimized. 

The maximum optical efficiency achieved reaches 88 % for a 3-QW structure. 

The width of optical intensity profile of MQW lasers increases, leading to lowering 
maximal light power density passing through laser facets, decreasing the risk of 
catastrophic damage of mirrors. 

It has been shown, by examining lasing offset voltage, that laser characteristic 
parameters are expected to be discontinuous functions of waveguide's width. 
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1.66 1.68 1.7 1.72 1.74 1.76 1.78 1.8 

U[V] 

Figure 6: Two Quantum Wells: Optical efficiency as a function of voltage applied for 
several combination of waveguides widths, as shown in the Figure: the first number is the 
width d w o, the second one is d w i (in nm). 



We did not study possible more complex structures. In particular, as modeling 
results show in case of 1-QW ([3] and [6j), introducing gradual doping profiles in 
waveguides, and also gradual changes to Al concentrations there may lead to signif- 
icant improvement of laser parameters. Also, the assumed here width of quantum 
wells, is close ony to the optimal one. 

Moreover, it is tempting to create structures with variable width of active re- 
gions (but the same lasing frequency). These would have, for instance, different 
characteristics as a function of temperature. 
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0.9 




U[V] 

Figure 7: Three Quantum Wells: Optical efficiency as a function of voltage applied for 
several combination of waveguides widths, as shown in the Figure: the first number is the 
width d w o, the second one is d w i (in nm). 
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Figure 9: Three Quantum Wells: Optical power intensity as a function of distance from 
the laser center, y. Calculations were done for the width d w \ = 0.130/xm, and several 
values of d w0 (in nm), as shown in the Figure. All curves are normalized to value of 1 at 
maxima. 
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